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ABSTRACT 


The economic feasibility of utilizing an ammonia bottoming cycle 
to improve the efficiency of a power plant is studied. The operating 
parameters of the bottoming cycle are examined to determine the optimum 
design binary cycle from a lifecycle costs viewpoint. 


The performance of the combined cycle is compared with that of a 
representative ‘modern’ steam plant: TVA's Bull Run. The analysis is 
based on current economic factors and the climatic conditions at the 
Bull Run site in Clinton, Tennessee. The optimum combined cycle is 
$98 million more profitable than the ‘Bull Run' plant over the expected 
lifetime of the plant. The combined cycle also consumes 2.7% less fuel 
than the unmodified steam plant. 


Optimization results are also presented for a range of environmental 
conditions at the plant location. 


The use of an ammonia bottoming cycle is economically as well as 


technically feasible provided the ambient water temperature is below 
69°F. 


gee 





ACKNOWLEDGMENTS 


I would like to take this opportunity to thank Professor 
A.D. Carmichael for his counsel and guidance which were indis- 
pensable during the course of this research. 

I also would like to thank my fellow graduate students 
in room 3-470 for their technical assistance as well as their 
friendship. 

I would like to extend my appreciation to the United States 
Navy for providing me with the opportunity to attend MIT under 
the auspices of the Junior Line Officer Advanced Educational 
(Burke) Program. 

Last, but certainly not least, I would like to express my 
gratitude to my wife, Lorraine for her patience and understanding 


during the research and particularly during the typing of this thesis. 


ae 





TABLE OF CONTENTS 


PAGE 
Chapter I INTRODUCTION iz 
Chapter II LITERATURE SURVEY 15 
Chapter III THE COMBINED CYCLE 17 
Chapter IV OPTIMIZATION PRELIMINARIES 28 
Chapter V THERMODYNAMIC ANALYSIS OF THE COMBINED CYCLE 50 
Chapter VI ECONOMIC ANALYSIS 62 
Chapter VII OPTIMIZATION PROCEDURE 6/ 
Chapter VIII OPTIMIZATION RESULTS 74 
Chapter IX SUMMARY AND CONCLUSIONS 92 
FOOTNOTES 95 
BIBLIOGRAPHY an 
Appendix I THE COMPUTER PROGRAM sk, 
Appendix I] STEAM CYCLE DATA 113 


whe 





Figure 
Figure 
Figure 
Figure 
Figure 
Figure 
Figure 
Figure 


Figure 
Figure 


Figure 


Figure 


Figure 
Figure 
Figure 
Figure 


Figure 
Figure 
Figure 
Figure 
Figure 


Figure 


CON DO FP WD NY — 


13 
14 
15 
16 


17 
18 


20 


2] 


22 


Figure 23 


LIST OF FIGURES 
Simplified Diagram of the Combined Cycle 
Heat Balance Diagram of TVA's Bull Run Plant 
Temperature-Entropy Diagram of TVA'’s Bull Run Plant 
Temperature-Entropy Diagram of Bottoming Cycle 
Temperature-Entropy Diagram of Combined Cycle 
Component and Flow Chart for Combined Cycle 
Temperature versus length for condenser 


Wetness Power Loss versus Steam Condensing Temperature 
for Steam Cycle 


Wetness Power Loss versus Steam Condensing Temperature 
for NH. Cycle 


Total Moisture Loss versus Steam Condensing Temperature 
for Combined Cycle 


Combined Cycle Efficiency versus Non-Dimensional 
temperature difference for various steam condensing temper- 
atures 


Temperature versus Length for Both zones of steam conden- 
ser/ammonia boiler 


Chen Factor versus Inverse of Martinelli-Nelson Factor 
Suppression Factor versus Two-Phase Reynolds Number 
Simplified Feedheater Flow Diagram 


Temperature-Entropy Diagram with Analysis Nomenclature 
for NH Cycle 


Optimization Procedure Flow Chart 


Condenser Water Velocity and Cycle efficiency versus 
Difference between NH. condensing temperature and water 
inlet temperature 


Combined cycle efficiency versus Non-Dimensional Tem- 
perature difference. Water inlet temperature = 55°F 


Steam Condenser/Ammonia Boiler Heat Transfer area versus 
Non-Dimensional Temperature difference 


Plant Lifetime Net Savings versus Non-Dimensional 
temperature difference. Water inlet temperature = 55°F 


Component Flow Chart for Combined Cycle 
Exit Losses Ammonia Turbine in Off-Design Conditions. 


abe 





Figure 
Figure 


Figure 
Figure 
Figure 
Figure 


Figure 


Figure 


24 
25 


26 


27 


28 


29 


30 


31 


Off-Design Performance of Combined Cycle 


Combined Cycle Efficiency versus Non-Dimensional 
temperature difference. Water inlet temperature = 35°F. 


Combined Cycle Net Savings versus Non-Dimensiona|l 
temperature difference. Water inlet temperature = 35"F 


Combined Cycle Efficiency versus Non-Dimensional 
temperature difference. Water inlet temperature = 45°F. 


Combined Cycle Net Savings versus Non-Dimensional 
temperature difference. Water inlet temperature = 45°F 


Combined Cycle Efficiency versus Non-Dimensional tem- 
perature difference. Water inlet temperature = 65°F 


Combined cycle Net Savings versus Non-Dimensional 
temperature difference. Water inlet temperature = 65°F 


Optimum Design Net Savings versus Condenser Inlet Tem- 
perature 





LIST OF TABLES 


Table 1 TVA's Bull Run 

Table 2 Description of Bottoming cycle T-S Diagram 

Table 3 Combined Cycle description 

Table 4 Summary of Engineering and Economic factors 

Table 5 Thermodynamic Specifications of the Combined 
Cycle 

Table 6 Optimum Design Parameters 

Table A2-1 Variation of Steam Cycle Efficiency and Tes with 


the Steam Condensing Temperature 





Ap 
AFLOW 


LMTD 


LIST OF SYMBOLS 
Area 
Boiling region heat transfer area 
Area perpendicular to flow 
Non-boiling region heat transfer for area 
Total heat transfer surface area 
Specific heat at constant pressure 
Specific heat at constant pressure 
Tube diameter 
Chen factor 
Friction factor 
Two-phase friction factor 
Acceleration of gravity 
32.2 ft-1bm/Ibf-sec® 
Enthalpy 
Heat transfer coefficient 
Modified enthalpy of vaporization 
Higher Heating Value 
Horsepower 
High-Pressure 
Interest rate 
Intermediate-Pressure 
Thermal conductivity 
Thermal conductivity 
Kilowatts 
Tube length 
Boiling region tube length 
Load factor 
Log mean temperature difference 


non 





LNB 
LP 


MW 


NH 
3 


PL 


p 
NH, 

PP 

psi 


PST 


Qn 
Qne 


Re 
Resp 
ReTp 


SC/AB 
TAV 
71) 


TF(2) 
5 ) 


Non-boiling region tube length 

Low Pressure 

Mass flow rate 

Megawatts 

Number of tubes in vertical row 
Number of years of plant life 
Ammonia 

Pressure 

Power loss due to turbine backpressure 
Ammonia cycle power 

Condenser Pumping power 

Pounds per square inch 

Steam cycle power 

Turbine power 

Heat transfer 

Heat added 

Heat transfer in boiling region 
Heat transfer in non-boiling region 
Heat rejected 

Reynolds number 

Single-phase Reynolds number 
Two-phase Reynolds number 

Nucleate boiling suppression factor 
Entropy 

Steam condenser/ammonia boiler 
Temperature 

Entropy averaged temperature of heat addition to ammonia 
Condensing temperature 

Ammonia boiling temperature 

Ammonia boiler inlet temperature 
Ammonia condensing temperature 


_9- 





Ah 
AP 
AT 
AT 
ATLM 
ATy 


"b 


Average temperature of heat addition 
Average temperature of heat rejection 
Ammonia temperature 

Steam condensing temperature 
Temperate-entropy 

Tube wall temperature 

Overall heat transfer coefficient 
Velocity 

Specific volume 

Pump work 

Thermodynamic quality 

Average quality in turbine stage 
Turbine exit quality 
Martinelli-Nelson factor 


SUBSCRIPTS 
Bull Run 
Combined cyle 
Extraction 
Property at saturated liquid state 
Property change associated with vaporization 
Property at saturated vapor state 
Property at saturated vapor state 


GREEK SYMBOLS 
Enthalpy difference 
Pressure difference 
Temperature difference 
Difference between Tc. AND TAV 
Log mean temperature difference 
Temperature difference between condensing fluid and tube wall section 
Auxiliary efficiency 
Boiler efficiency 


ee 





Cycle efficiency 

Generator efficiency 

Turbine isentropic efficiency 
Mechanical efficiency 

Overall cycle efficiency 

Pump efficiency 

Steam cycle thermal efficiency 
Steam cycle thermal efficiency 
Thermal efficiency 

Dynamic Viscosity 

Average dynamic viscosity 
Density 

Surface tension 


== 








IT. INTRODUCTION 

During the past few years, as the demand for electrical power has 
increased, mankind has become increasingly concerned over several 
important factors which will have impact on how electrical power will 
be generated in the future. The first of these considerations is the 
dramatic increase in the price of fuel and the realization that our 
fuel resources are limited. Secondly, we have become more concerned 
over how the electrical power generating plant affects our environment. 
The specific environmental considerations are how the fuel combustion 
process pollutes the air and how the heat which the plant rejects affects 
the temperature of the heat sink, usually a lake or river. 

In order to partially overcome these problems a plant with a higher 
thermal efficiency than presently exists is desired. The higher effic- 
lency plant would, for the same power output, burn less fuel and there- 
fore exhaust fewer combustion products to the atmosphere. Furthermore, 
as the plant efficiency increases the amount of heat rejected per unit 
power output decreases. 

The problem is: How to increase the efficiency of generating 
electric power? 

There are many ways to design a highly efficient power plant. Some 
of the proposed 'advanced' cycles are both expensive and exotic. The 
possibilities include: 

(1) Fuel-cells 

(2) MHD 

(3) Combined steam/gas turbine cycles 
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(4) Liquid metal topping cycles. 
(5) Bottoming cycles utilizing a refrigerant-type working fluid. 
In this thesis, the bottoming cycle will be examined. 


The thermal efficiency of an ideal cycle n,, = 1-7 (Equation 1) 


where Ty is the average temperature at which heat is rejected and Ty, 
is the average temperature of heat addition to the cycle. 

The purpose of a bottoming cycle is to increase the thermal efficiency 
of a cycle by lowering Tj. Conventional steam plants are unable to take 
full advantage of a low heat rejection or condensing temperature due to: 


(1) The large specific volume of steam at low condensing 
temperatures and pressures. 


(2) The increased moisture content in the form of droplets 
as the condensing temperature drops. 


The above considerations require large, costly, and relatively 
inefficient low pressure turbine stages. 

In order to overcome these limitations and to take full advantage 
of a low temperature heat sink, another Rankine cycle using a different 
working fluid could be added to a steam plant. This 'bottoming' cycle 
would receive the energy which the steam plant rejects and would have 
the capability of lowering condensing temperatures. 

The selection of a working fluid for the ‘bottoming’ or sub-position 
cycle cannot be accomplished by a closed form equation, but rather by a 
series of trade-offs. The desired properties of the working fluid are: 


(1) Have a relatively low specific volume at the condensing 
temperatures to be considered. 


(2) Be a liquid at ambient temperatures, so it can be pumped. 
(3) Be non-toxic and non-corrosive. 
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(4) Have a low specific heat and a large latent heat of 
vaporization over the temperature range to be con- 
sidered. (This results in the average temperature of 
heat addition being as high as possible). 

(5) Have a condensing pressure above atmospheric pressure. 

(6) Have a high speed of sound at turbine exit conditions. 

(7) Have good heat transfer properties. 

(8) Be chemically stable and non-flammable. 

(9) Be inexpensive and plentiful. 

No known substance satisfied all of these requirements, but 
Ammonia (NH3) seems to be the best overall choice. 

The feasibility of the ammonia bottoming cycle will be evaluated 
by investigating a representative ‘modern’ steam plant modified to 
include the sub-position cycle. The thermal efficiency, capital, 
and operating costs of the binary cycle plant will be compared with 
those of the unmodified 'modern' plant. With these factors in mind 
an optimum design combined cycle plant will be generated, which will 


result in maximum savings to the owner/operator over the expected 


lifetime of the plant. 
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II. LITERATURE SURVEY 

The concept of a bottoming cycle is not new. In 1961, Aronson” 
proposed such a binary-vapor cycle, utilizing Freon -12 as the working 
fluid. His analysis showed a heat rate improvement of 1-2% over that 
of a 'modern' steam plant as well as establishing the economic feasi- 
bility of the cycle. 

In 1969, Wood? stated that a binary-vapor bottoming cycle was not 
feasible, citing available energy loss due the temperature difference 
across the steam condenser/secondary fluid boiler. 

The Aronson proposal has been supported by Steigelman et al’, who 
demonstrated the desirability of a binary-cycle using cooling towers 
provided the ambient air temperature was low enough. Furthermore, £1 
Ramly and Budenholzer” professed the advantages of a bottoming cycle to 
Operate in conjunction with a nuclear power plant. 

Perhaps the most comprehensive study to date has been that of 
Slusarek®, in which the feasibility of the ammonia bottoming cycle 
was studied in detail. The Slusarek report included a thorough component 
design, an economic analysis, and an engineering optimization. While 
the majority of the report is accurate, two significant shortcomings 
must be noted: 


(1) The heat transfer coefficients in the ammonia boiler were 
calculated using inappropriate correlations. 


(2) The impact of the steam condensing temperature and the 
temperature difference between the main and bottoming cycles 
were not adequately examined. 
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In this study, the operating parameters of a bottoming cycle will 
be studied and optimized in order to prepare a preliminary design of the 


best plant from an engineering as well as economic viewpoint. 
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III. THE COMBINED CYCLE 

As outlined in chapter one, the bottoming cycle utilizes the heat 
rejected by a steam plant as its energy source. A simplified diagram of 
a combined steam-ammonia plant, hereafter called the ‘combined’ cycle, 
is shown in figure one. 

To carry out the analysis of the feasibility and economic optim- 
ization of the combined cycle, an existing ‘modern’ steam plant was modified 
to include a bottoming cycle. The Tennessee Valley Authority's ‘Bull 
Run' plant was selected for this role. 

Bull Run is a 900 MW, coal fired, supercritical steam plant which 
first went into operation in the mid 1960's. The plant features a 
boiler exit temperature of 1000°F and a single reheat to the same tem- 
perature. The design condenser pressure is 1.5" of mercury, which 
corresponds to 91.7°F and .74 psi. The thermal efficiency of Bull 
Run, excluding the boiler losses is 45.78%. Including the boiler 
efficiency of 89%, the cycle efficiency is 40.74%.7 

The Bull Run plant flow diagram and the Temperature vs. Entropy 
diagram are shown as figures two and three respectively. 

The ammonia bottoming cycle is a simple Rankine cycle. Feedheating 
must be utilized since the feasibility analysis is sensitive to small changes 
in the efficiency of the bottoming cycle.® 

Based on earlier studies, a four stage turbine for the ammonia cycle 
seems to be a logical choice. Consequently, the cycle will have three stages 
of feedheating. Data from steam plants indicates that the maximum efficiency 
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TEMPERATURE (2 ) 


ENTROPY 


FIGURE 3 TEMPERATURE - ENTROPY DIAGRAM OF TVA'S BULL RUN PLANT 
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increase for three stage feedheating occurs when the working fluid 
temperature is elevated three-quarters of the difference between the 
condensing and evaporating temperature by feedheating.~ 

The Temperature-Entropy diagram of the bottoming cycle is shown 
as figure 4. Note that the ammonia is not superheated. In order to 
maximize the efficiency of the combined cycle, the average temperature 
of heat addition to the ammonia cycle must be as close as practical to 
the steam condensing temperature. If Superheat were introduced, the 
average temperature of heat addition for the bottoming cycle would 
drop, negating the efficiency gains usually derived from superheating. 

In order to bring the steam and ammonia cycles together into one 
combined cycle certain modifications must be made to the two cycles. 
The modifications are: 


(1) Replacing the water-cooled steam condenser with a steam 
condenser/ammonia boiler. 


(2) Raising the backpressure of the low pressure steam turbine 
and combining the remaining LP turbine stages with the IP 
turbine. 


(3) Removal of feedheating stages in the steam plant made un- 
necessary by raising the steam condensing temperature. 


The T-S diagram and component flow chart for the combined cycle 
are shown as figures five and six respectively. 
Since the high temperature end of the steam plant is not altered, 
the only independent parameter in the steam side of the combined plant 
is the steam condensing temperature, To. On the ammonia side there are two 
independent variables: 
(1) The ammonia boiling temperature TF(1). 


(2) The ammonia condensing temperature, TF(5). 
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FIGURE 4 T-S DIAGRAM OF BOTTOMING CYCLE 


DESCRIPTION 

Condenser exit (saturated liquid) 

Exit of third feedheater, boiler inlet (liquid) 
Boiler exit, turbine inlet (saturated vapor) 
Turbine exit, condenser inlet (liquid and vapor) 
DESCRIPTION 

Mass flow rate through boiler 

Mass flow rate through condenser 

Mass flow rate of feedheating extraction #1 
Mass flow rate of feedheating extraction #2 
Mass flow rate of feedheating extraction #3 


TABLE 2 DESCRIPTION OF T-S DIAGRAM (FIGURE 4) 
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FIGURE 5 T-S DIAGRAM OF COMBINED CYCLE 
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TABLE 3: 


COMBINED CYCLE DESCRIPTION 


STREAM NUMBERS REFER TO FIGURES 5 AND 6 


STREAM 

Fuel 

Air 

Exhaust 
Water 
Supercritical Liquid 
Steam 

Steam 

Steam 
Steam/Water 
Ammonia 
Ammonia 
Ammonia 
Ammonia 
Water 


Water 


* INCLUDES FEEDHEATING 


NO. DESCRIPTION 

a Fuel inlet 

b Air inlet for combustion 

C Combustion products 

] Condenser exit 

Ze Boiler inlet 

3 Boiler exit/HP turbine inlet 

4 HP turbine exit/reheater inlet 
° Reheater exit/IP turbine inlet 
6 IP turbine exit/condenser inlet 
7 Condenser exit 

8* Boiler inlet 

9 Boiler exit/NH3 turbine inlet 
10 NHz turbine exit/condenser inlet 
1] Condenser coolant water inlet 
4 Condenser coolant water exit 
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The three variables mentioned above are the only independent 


variables in the analysis. The significant dependent variables are: 


(1) The entropy averaged temperature of heat addition to the 


(2) 


ammonia cycle, TAV, which is a function of TF(1) and TF(5). 


The difference between the steam condensing temperature T., 
and TAV. This variable is called AT and is a function of 
all three independent variables. The magnitude of this 
parameter is related to the available energy lost due to 
the temperature drop between the cycles. 


The Log Mean Temperature Difference (LMTD) in the steam 
condenser/ammonia boiler. This factor which is numerically 
close to AT, is a function of the same three variables, 

and is a significant consideration in the sizing and cost 
estimation of the steam condenser/ammonia boiler. 


The object of the combined cycle optimization is to consider 


how the three independent variables affect the plant performance and 


cost. 
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IV. OPTIMIZATION PRELIMINARIES 

In this chapter the impact of the independent and dependent 
variables on the analysis will be examined. In addition, limita- 
tions on the range of these parameters will be discussed where 
appropriate. 

a. The ammonia condenser and the ammonia condensing temperature. 

In order to maximize the overall cycle efficiency, the lowest 
possible ammonia condensing temperature, TF(5), is desired. However, 
for a power plant with a water-cooled condenser there are three factors 
which limit how low the condensing temperature can be. The factors 
are: 


(1) The condenser coolant water temperature at inlet must be 
greater than 32°F. 


(2) The water velocity through the condenser must not be so 
high as to cause excessive corrosion in the tubes. 


(3) The amount of power required to pump the coolant water 
through the condenser tubes must not negate the advantages 
of lowering the condensing temperature. 
The reason why the water temperature at inlet must be greater 
than 32°F is that it must be a liquid in order to be pumped. If the 
ambient temperature is below 32°F, a cooling tower would be utilized 
in place of a wet condenser. The considerations of water velocity and 
pumping power are not quite so obvious, but are related to each other. 
In a heat exchanger, such as a condenser, the temperature difference 
which will partially determine the condenser size is 
nie= 1. - (Ty + Tg)/2 (equation 2) 
where Tr is the condensing temperature, T; is the water inlet temperature, 
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and Tp is the water exit temperature. The temperature vs. length plot 
for a condenser is shown figure 7. 

The heat transferred in the condenser Q may be expressed as: 
Q = UA AT (equation 3) 
where A is the heat transfer area and U is the overall heat transfer 
coefficient for the heat exchanger. The overall heat transfer coefficient 
may be computed by combining the thermal resistances in the unit as if 
they were resistors in a parallel electrical circuit. Assuming the 
thermal resistance of the condenser tubes to be negligible, and the 
thermal resistances for convection and condensation to be the 


inverse of the average heat transfer coefficients on their respective 


Sides of the unit, then U = 1/(1/Ncondensing * '/Mconvection on) (equation 4) 
fluid water side 
For shorthand purposes the following symbolism will be used: 


Ncondensing fluid > NNH3 cond and hconvection on water side > Mwater 
In terms of the fluid and thermal parameters for a horizontal tube 


condenser: 


NH = 728 [1+.2C AT(n-1)/hfg] “Ag(o)(o - py)k3 h'¢g)/(nDu ATy) (equation 5) 
con 


Nwater = -023 k¢/D (p¢VD/u¢) °8 (Cpu/kr) °4 (equation 6) 
where [ = specific heat of condensing liquid 

Cy = specific heat of water 

D = tube diameter 

g = acceleration of gravity 

hég = heat of vaporization of condensing fluid at condensing 

pressure 

h' ¢g = 1+(.68 C IN Neg 

k = thermal conductivity of condensing liquid 

Ke = thermal conductivity of water 
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n = number of tubes in vertical row 
le) = density of condensing liquid 


Pe = density of water 


P, = density of condensing vapor 

u = dynanic viscosity of condensing liquid 

Us = dynamic viscosity of water 

V = water velocity through the tubes 

AT, = temperature difference between condensing fluid and 


tube wall 
Furthermore, an energy balance for the condenser states: 
Q=m water Oat a (Ty 7 1) 7 "condensing fluid *e Neg (equation 7) 
and solving for To ~ Ty» the water temperature increase: 
To - Ty =(m condensing fluid *e heg)/ ("water eDeetey (equation 8) 
where X, is the condensing fluid quality. Combining equations (2), (7), 
and (8) then, 
ee’ - (1; + Tooeiea= IC = iietO/ citar an) (equation 9) 
From equation (3) it may be seen that for a condenser to be of 
minimum size and cost once the heat transfer amount Q has been selected, 
the product UAT must be as great as possible. 
However, from equations (5) and (6), it may be shown that AT eances 
much more rapidly than U, apparently stating that AT sholildsbe maximized. 
Equation (9) indicates that once Tc is chosen (and ideally as low as possible), 
there are only two ways to increase AT. The first method is to decrease 
Ty, the water inlet temperature, which is difficult since Ty; is a function 
of the plant site. 
On the other hand, if the Q/2myatey term in equation (9) is small 


the value of AT will remain large. The problem is that the only way to 
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decrease the Q/2m, ater term for a fixed Q, is to increase the mass 

flow rate of water through the condenser tubes. This brings one 

back to the problems of water velocity and pumping power stated earlier 
in this chapter. 

In order to keep the economic analysis of the components as 
Simple as possible, the combined cycle plant condenser under consider- 
ation will have the same ratio of cycle power output to surface area 
as the 'Bull Run' unit. This assumption is valid, since the condenser 
limitations are on the water side. The ‘Bull Run' condensing system 
consists of 41368 tubes, each of which has an outer diameter of .875", 

2 


and a total heat transfer surface area of 320,000 feet™. Since the 


mass flow rate, Mater = Pwater Vwater Aflow (equation 10) 


and A¢jow = nnd°/4 (equation 11) 
(where Ag], is the area perpendicular to the direction of water flow, 


n is the number of tubes, and D is the tube inner diameter), then 


Vwater = Twater/Pwater Aflow = 4twater/owater nD” (equation 12) 
Now that an expression for the water velocity as a function of 
the mass flow rate and the condenser geometry has been defined, a limit 
on the mass flow rate and TF(5) can be specified once a limitation on 
the water velocity is established. 
Since ammonia reacts readily with copper, stainless steel is the 
likely candidate for the tube material. Studies indicate that for steel 
tubes corrosion will not be excessive as long as the water velocity is 


below 15 oy 
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This water velocity may require excessive pumping power which 
would be counter-productive. The pumping power, PP may be expressed 
as: PP = (AP x Q)/n, (equation 13) 
where Q is the volumetric flow rate, Np the pump efficiency, and 
AP is the pressure drop through the tubes. 

The volumetric flow rate Q, is 
Q = ™ water/? water (equation 14) 
and the pressure drop, AP is: 
fee 4 £(L/D)(oyater Vwater/290) (equation 15) 
where f is the pipe friction factor, L the length and Do is the 
units correction factor. 

For turbulent flow in tubes: 
fo A079/R, °° (equation 16) 
where R, is the Reynolds number and is equa! to 


R (equation 17) 


e- (Pater Vwater 2)/ Ywater 
Combining equations (13) through (17), the total pumping power 


1S: 


2 : 
PP = .316 / (o,V,D/u,)-29 (L/D) (© V 72g.) (mh /oy)(1/ny) (eq. 18) 


recalling that L = TOTAL HEAT TRANSFER AREA (A)/(nmD) (eq. 19) 
and m= (OV nnD*)/4 (equation 20) 
the pumping power expression simplifies to: 

Bet.079 Au” p°!? vo-!y7 (4 0°? g, np) (equation 21) 

in which V is the most significant factor and the only one that will 


change noticeably from one set of parameters to another. 
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Thus far it has been shown how the ambient water temperature 
re and practical engineering considerations limit the ammonia 
condensing temperature. 

To summarize: 


(1) Once a plant site has been chosen 
VN is fixed. 

(2) For a fixed T; and a given condenser, 
the mass flow rate (and velocity) must 


be maximized in order to keep the ammonia 
condensing temperature as low as possible. 


(3) Limitations on the water velocity are 
imposed by corrosion considerations. 


(4) In addition to #3, the water velocity must 
not be so large as to require excessive pumping 
power which would negate the advantages of 
lowering TF(5), the NH, condensing temperature. 
The interface between the gains of lowering TF(5) and the 
pumping power requirements will be considered later. 
b. Steam Condensing Temperature, T, 
It is desired to find the steam condensing temperature which 
will give the maximim cycle efficiency. To obtain the optimum To, 
it must be determined whether the steam or ammonia cycle makes better 
use of the low temperature end of the cycle. 
Assuming that the steam and ammonia turbines have the same 
isentropic efficiencies in the two-phase region, the deciding factor 
in the selection of To, is the power loss due to moisture droplets in 


the turbine. 


Common practice indicates that the power loss due to wetness is 
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one percent for each percent of average wetness in that turbine stage. 

Applying this wetness consideration to the ‘Bull Run' plant, a 
plot of power loss due to wetness versus the steam condensing temperature 
may be created as shown in figure 8. 

Once the ammonia condensing temperature TF(5) has been selected, 
the power loss analysis may be carried out for the ammonia cycle as 
well. However, since the third independent variable, the NH, boiling 
temperature TF(1), has not been fixed the wetness calculation must be 
carried out for each value of TF(1). This procedure produces a family 
of curves one of which is shown in figure 9. 

Figure 8 demonstrates that the power loss for the steam cycle 
increases aS Ts drops, and figure 9 shows the opposite trend for the 
ammonia cycle. 

The optimum value of Ty may be found by superimposing figures 8 
and 9 and plotting the combined cycle wetness loss versus Ts. 

Figure 10 shows the superimposition of figures 8 and 9 and the 
resulting total wetness loss. Note that a mimimum loss occurs around 
ic = 131 F for the parameters chosen. Furthermore, if the combined 
cycle efficiency is plottedversus the non-dimensionalized temperature 
difference between the steam and ammonia cycle for various values of 
Tg; as shown in figure 11, it is clear that To= 131 F is the optimum 
choice. 

The selection of a steam condensing temperature at this level is 


different from the studies cited in chapter II. All of the earlier work 
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assumed a value of Tc in the neighborhood of 220° - 250°F without 
considering the wetness losses. 

The analysis which generated figures 10 and 11 will be discussed 
in the next chapter. 


c. The Ammonia Boiling Temperature TF(1) and the Steam Condenser/ 
Ammonia Boiler. 


The third independent variable is TF(1), the ammonia boiling tem- 
perature. This variable is most significant in its impact on the cycle 
efficiency and the size and cost of the steam condenser/ammonia boiler. 
The value of TF(1) is not as important as its relationship with the steam 
condensing temperature Tc. This relationship defines the dependent 
variable AT, the average temperature difference between the steam and 
ammonia cycles. 

As indicated in chapter three and graphically shown in figure Il, 
the combined cycle efficiency is highest when AT js as small as possible. 
Unfortunately, since the average temperature difference AT is numerically 
Similar to the log mean temperature difference in the steam condenser/ 
ammonia boiler and the heat transferred in the heat exchanger Q is: 


Q = UA AT (equation 22) 


Log mean’ 
Thus as the temperature difference decreases the heat exchanger area 
increases. 

It is obvious that the efficiency gains derived by a small AT may 
be negated by the additional capital cost of a large steam condenser/ 


ammonia boiler. 
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The task at hand is twofold: 
(1) Determine how AT affects the overall thermal efficiency. 


(2) Determine how AT affects the steam condenser/ammonia 
boiler size and cost. 


Task one will be considered in the next chapter. 

In order to determine how AT affects the heat exchanger size and 
GOSt, a detailed heat transfer analysis of the unit must be made. 

The steam condenser/ammonia boiler to be considered is a shell 
and tube arrangement with the steam condensing outside the array of 
horizontal tubes through which the ammonia flows. 

For analytical purposes the SC/AB may be divided into two zones. 
In the first zone, the ammonia is heated from the feedheating outlet 
temperature, TF(2) to the ammonia boiling temperature TF(1). The heat 
transfer mechanism on the ammonia side in this region is forced con- 
vection, with subcooled boiling effects being neglected, and on the 
Steam side the mechanism is condensation. 

In the second zone, the ammonia temperature remains constant as 
the ammonia is heated from a saturated liquid to a saturated vapor state. 
The heat transfer mechanism on the NH side in this zone is very complex, 
and includes effects of nucleate boiling and forced convection. The heat 
transfer mechanism on the steam side is once again condensation. 

In both the non-boiling (zone 1) and boiling (zone 2) regions the 
heat transfer coefficient on the steam side Noteam May be expressed as: 


= .728 [1+.2C AT(n-1)/h¢g] Ag(6) (0 Py)k h'¢g)/nDu AT) @q. 23) 


where the equation and symbols are identical to equation 5, which was used 


h 
steam 


Shale 





ZONE II 


TF(1 
TEMPERATURE 


AMMONIA 


TF(2 


NON- BOILING 
ZONE BOILING ZONE 





Ie— Lng —>————-_ | ——-——_-—> 


LENGTH 


FIGURE 12 TEMPERATURE VERSUS LENGTH FOR BOTH 
ZONES OF THE STEAM CONDENSER/AMMONIA 
BOILER 


Sees 








in the NH. condenser calculations. 


On the ammonia side in the non-boiling region the forced convection 


i . ae, 4 
coefficient is: hyy, p, = .023 (k¢/D) (o¢VD/u) (Cou/k¢) (eq. 24) 


which is the same expression as equation(6),except that the fluid 
properties for ammonia are used here. 

In zone 2, the boiling region, the heat transfer coefficient may 
be calculated by using several empirical relationships. The most 
reliable seems to be the Chen correlation |¢ which states that the 
heat transfer coefficient in the two-phase region where boiling and 
forced convection are present, may be calculated as follows: 

HTWO PHASE = "NUCLEATE BOILING X S + heoRCED CONVECTION 2 PHASE XF (eq. 25) 
where S is the boiling suppression factor, a measure of how much 

the fluid motion upsets the normal boiling mechanism and F is the 

Chen factor which is a measure of how much the fluid motion induced by 
boiling upsets the normal forced convection mechanism. 

The Chen factor may be calculated as follows: 


(1) Determine the Martinelli-Nelson factor X++, where 


: 5 j 9 

Xet (24/0 ¢) (U¢/u,) ((1-x)/x) (equation 26) 
Note that po, and op, represent the density of saturated vapor 
and liquid respectively, U¢ and u, the viscosity of saturated 


liquid and vapor respectively, and x is the thermodynamic 
quality. 


(2) Utilizing figure 13, which displays F versus /X, read off F. 
The boiling suppression factor S may be found as follows: 


(1) Find the two-phase Reynolds number Rel which is represented by 
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(3) Utilizing steps one and two the forced convection of ammonia 
heat transfer coefficient may be found. 


(4) Solve the heat flux equation by 
QnB/Ang = steam (Ts-Twa1l) = MNH3 Fo (walt-TnH3) (eq. 33) 
iteration. The reason that this relationship cannot be solved 


in closed form is that the an for heteam (equation 23) 


contains the term (Ts-Tya11) “making a numerical result 


preferable. The iteration is carried out by assuming a tube 


wall temperature and incrementally changing 1,51] until Deena 


(Ts-Twal1) equals ys (Twa11-"NHg)- Note that Tyy3 is the 
numerical average between the ammonia inlet and outlet temperature. 


(5) Divide the non-boiling heat transfer by the heat flux to 
determine the area of the non-boiling section, Aye. 


(6) Calculate the pressure drop in this section of the SC/AB and 
determine that the design choices made in step one are 
responsible. Use equation 15 and substitute ammonia properties 
for those of water. 

The procedure for the boiling section calculations are basically the 

Same, but are more complex. 

Note that in equation(26)and(30),the variable 'x' which represents 
the thermodynamic quality appears. In order to calculate the heat 
transfer coefficient for a region where the quality goes from zero to one, 
what value of x should be used? A simple average between the end values 
is not suitable since in both equations(26)and(30),the 'x' term is non- 
linear. However, in order to make this problem less difficult to solve, 
an approximation which yields good accuracy is to use an average quality 
in the analysis, provided the change of the quality is small. 


In keeping with this idea, the boiling region of the SC/AB is 


divided into ten segments, each having a quality change of one-tenth and 


ye 





and utilizing the average quality in that segment in equations (26) and 
30) . 

The configuration of the steam condenser/ammonia boiler is determined 
(step one for non-boiling section) and the amount of heat transfer in 
the boiling region, Q, » as well as the NH, mass flow rate (step 2) 
are fixed. Then for each quality segment of the boiling region the heat 


flux iteration must be performed: 


Q : 
B SEGMENT = h 43 (To-T 


B SEGMENT wa MBSE TREHESEcony PFO) Tarn) (ea. 98) 
This iteration is more complex than the non-boiling case since hynycLEATE 
BOILING 
and eam both contain temperature difference terms raised to fractional 
powers. 


Once the heat flux for each segment has been determined, the heat 
transfer area of that segment may be found. The total boiling region area 
Ap is the total of the areas of the ten quality segments. 

Finally, the pressure drop in the boiling region must be evaluated 


to determine design feasibility. Using a homogeneous flow Welle. the 


pressure drop for the boiling region is: 


2 ; 
AP = 2f+p pVL 1+ Veg 4 ov" (Vfq/vF) (equation 35) 
D 2Ve 


where frp = .079/ (sly oo 


(equation 36) 
and u is the average viscosity and may be expressed by the McAdams cor- 


Peeeion |’: T = Xug + (1-x)u.. (equation 37) 
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Once again quality is important and the calculation of the pressure 
drop is carried out in the same manner as the heat transfer coefficients, 
by assuming the calculation to be piecewise linear. 

d. Chapter summary 

In this chapter, the impact of the three independent variables 


TF(5), TF(1), and To on the heat exhanger design and cycle efficiency 


have been examined. Specifically, the NH, condenser performance and 
pumping requirements, the selection of an optimum To, and the sizing 

of the steam condenser/ammonia boiler have been considered. In the next 
chapter, the thermodynamic analysis will be discussed and the impact of the 


three variables on cycle efficiency. 
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V¥. THERMODYNAMIC ANALYSIS OF THE COMBINED CYCLE 

In this chapter, the impact of the three design parameters, the 
ammonia boiling temperature, ammonia condensing temperature, and steam 
condensing temperatures on the combined cycle efficiency will be studied. 

a. The Steam Cycle 

The steam condensing temperature, Tos is the only variable which 
affects the thermal efficiency of the steam cycle. Specifically, it 
would be desirable to know how the steam cycle efficiency Not» varies 
with Tos the steam condensing temperature. 

The power output of the steam cycle is affected by two considerations 
as the steam condensing temperature is increased to accomodate the bottoming 
cycle. The first of these considerations is that the LP steam turbine 
produces less power due to the higher exhaust backpressure. Secondly, 
since the steam is being condensed at a higher temperature than the 
‘Bull Run' plant, less steam must be extracted from the power producing 
turbines for feedheating purposes. Note that those considerations 
affect the cycle efficiency in opposing manners. 

With regard to the loss of power due to increased backpressure, 
only the low pressure steam turbine is affected. The power loss associated 


with the higher backpressure is: PL = Motoam ( hi - ho ) (equation 38) 


where h; is the enthalpy of the stage exhaust steam for the increased back- 
pressure case and hy is the stage exhaust enthalpy for the unmodified 


Bull Run plant. Equation (38) must be modified to include the difference 


-50- 





in wetness losses between the original and the modified steam plant. 
Accordingly if. h. is the inlet enthalpy, x, is the average quality for © 
the modified plant, and xX, is the average quality for the base case, | 


equation (38) becomes: 
PL = mM steam ((hy - No) %> - (hy ~ 9) %) (equation 39) ey 

Equation (39) is valid for each stage in the LP turbine provided . 
that the proper steam flow rate, M oteam is used. 

In addition to the difference in power output due to the higher 
backpressure, the changes in the feedheating extractions also affect the 
power output. 

In the ‘Bull Run’ plant, the extracted steam is condensed in a 
closed-type feedheater, transferring heat to the feedwater until the ® 
condensed extraction temperature is 5°F higher than the feedwater outlet 


temperature. Ultimately, the condensed extraction is mixed with the Pp 
feedwater from the main steam condenser. ‘i 
In order to determine the cycle power variation due to changes in 
feedheating extractions, an expression for the mass flow rate of steam 
‘a’ in figure 15 must be established. 
Referring to figure 15, an energy balance of the closed feedheater 


indicates: 


Ma (ha - hy) = ty (he - hg) (equation 40) , 
A design constraint for the ‘Bull Run' feedheaters states: 


meet, + o°F (equation 41) 


i 
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and since C = 1 BTU/Ibm - °R, then: 


pwater 
hh = ha + 5 BTU/1bm (equation 42) 
An energy and mass balance in the mixer states: 
ig = ma + te (equation 43) 
hg = (tig hy + moh. )/ (ing + me) (equation 44) 
Combining equations (40) through (44), the extraction flow, m,, 
may be expressed as: 
my = M.((ing + M)he - Mg(hy + 5) -thche 
(wh, + me)(ha - 2he - 5) + ma(hg + 5) + moh, |(equation 45) 
The condenser flow rate is also affected by the amount of 


feedheating extraction. The definitive relationship is: 


M. = MBOTLER - "QTHER EXTRACTIONS ~ ™a (equation 46) 
renaming MgoILER - "OTHER EXTRACTIONS = "Wg (equation 47) 
and combining equation (45), (46), and (47): 

m, = (m, - m, ) m 


a By SUNY (m, aR 


g (eq. 48) 
m (ha - 2h, - 5) + m(h, +5) + (m, - m_ dh 

The variables in equation (48), M, and i are fixed at the values 
of the unmodified plant. 

The feedheater exit enthalpy will be fixed as a design constraint 
in this analysis. Thus, there are only two unspecified variables, the 
first is the enthalpy of the main condenser condensate, h_» which is a 
direct function of the steam condensing temperature, Ty. The second 
variable is the extraction mass flow rate, m,. Based on the arguments 


presented, it is now possible to state that: my = f(T.) (eq. 49) 
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where the function f is specified in equation (48) and in the 
relationship between h. and To» which may be found in the steam 
tables. 

Therefore, if the extraction mass flow rate for the unmodified 
cycle is Ma and the extraction flow rate for the new steam con- 
densing temperature is Ma, » the power loss due to the higher con- 
densing temperature (including the effects of equation (39)) is: 
PL = Ms Ch, - ho)X, - me (hy : hx, (equation 50) 
where m, is the turbine stage mass flow rate for the unmodified 
case, and mM, is the modified cycle mass flow rate. The modified 
cycle flow rate may be expressed as: 


MS] a SO "4 ue 


Combining equations (50) and (51), the power loss is : 


: m. ) (equation 51) 


PL = mg Ch, = ho) Xo - (Ms, x (ma | = Ma, ) (hj ~ hy) oe veg: 62) 
The thermal efficiency of the modified steam cycle becomes: 


Nth BCP = PLY Qn eam (equation 53) 


Steam 
where P, is the unmodified cycle power output and Ure is the heat 


added to the steam cycle. 


Consequently, the heat rejected by the steam cycle, which is 


also the heat added to the bottoming cycle is: 


QR = Qa, = QA (l-n 7 
equation 54 
steam NH3 steam Chives) (eq ) 
It may be noted that effect of only one stage of feedheating was 
shown in the modification of the 'Bull Run' plant. For the range of 


values where the steam condensing temperature Ts yields optimum combined 
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cycle efficiency, only the first feedheating stage is affected. 
However, a more complex analysis was undertaken for values of T, 
which affected the second and third feedheating stages. 

b. The Ammonia Bottoming cycle 

The ammonia bottoming cycle is a Rankine cycle with three stages 
of feedheating and no superheating. 

In carrying out the thermodynamic analysis of the bottoming cycle, 
the following symbols are used: 


S = entropy (BTU/1bm °R) 


h = enthalpy (BTU/1bm) 


p 


Pressure (psi) 


TF = temperature (°F) 


v = specific volume (ft?/1bm) 
f = used as a subscript, indicates saturated liquid state 
g = uSed as a Subscript, indicates saturated vapor state 


fg = used as a Subscript, indicates property change associated with 
vaporization or condensation 


Figure 16 is the T-S diagram for the bottoming cycle with the no- 
menclature to be used in the analysis. Note that since there are 
three stages of feedheating, there are five temperature levels in the cycle. 
The temperature levels are numbered sequentially, starting with level one 
which is the ammonia boiling temperature and culminating with level five 
which is the ammonia condensing temperature. 

The power extracted from each turbine stage is the product of the 
mass flow rate and the enthalpy drop in the stage, minus mechanical, kinetic 


energy and moisture losses. 
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The ideal enthalpy drop is found by considering the expansion 
in the turbine stage to be isentropic. Therefore, using the first 
stage as an example: 
s(2s) = s(1) (equation 55) 
and the quality at state 2s is: 
x(2s) =(s(2s) - $¢(2))/S¢q(2) (equation 56) 
therefore, the enthalpy of state 2s is: 
h(2s) = h¢(2) + h¢g(2) mx( 2s) (equation 57) 
However, since in reality the expansion is not isentropic, nj., the 
isentropic efficiency is introduced: 
Nig =(h(1) - h(2))/(h(1) - h(2s)) (equation 58) 
therefore the actual enthalpy drop in the first stage is: 
Ahjo = h(1) - h(2) = Mean(l} = h(2s)) (equation 59) 

The loss due to wetness corresponds to one percent power loss 
for each one percent of average moisture in that stage. Conversely, 
the wetness losses could be accounted for by multiplying the enthalpy 
drop by the average quality in the stage. For the first stage,the exit 


quality x(2) is: 


eee=(h{2) - he(2))/ hg(2) (equation 60) 
and the average quality in the first stage 
xy = (x(1) + x(2))/2 (equation 61) 


The mechanical, stage, and exit losses are each assumed to be one 
percent. Expressing these three factors in one expression called 
mechanical efficiency, n, which is simply one minus the three one 


percent loss, thus its =, Clee (equation 62) 
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Therefore, the power extracted from the first stage of the ammonia 
turbine PT? is: 

Pr, = my(h(1) - h(2)) - ny > (x01) + x(2))/2) (equation 63) 

The same analysis may be applied to the final three stages of the 
ammonia turbine, provided the proper mass flow rate for each stage is 
used. Recall that each successive stage has a smaller flow rate due 
to feedheating extractions. 

The purpose of feedheating is to improve the cycle efficiency by 
increasing the average temperature of heat addition to the cycle. 

This is accomplished by extracting some of the working fluid and 
utilizing it to heat the boiler feed. The type of feedheater used in the 
analysis iS an open-type, in which the feed is mixed with the extracted 
fluid. 

The extracted fluid mixes with the feed after the feed pump exit. 

The mass flow rate of the extraction is adjusted so that the mixture is 

at the saturated liquid stage. For example, the condensate (state F5 of 
figure 16) is pumped to the pressure corresponding to saturation at tem- 
perature level four. Ammonia is extracted at state 4 and mixed with the 
feed until the mixture is at the state F4. This process is repreated until 
the feed is at the boiler inlet condition. 


An energy balance in each feedheater states: 


MFeed AMteeq = Mextraction Shextraction (equation 64) 
or Mextraction ~ "feed Ahfeed/ Ahoytraction (equation 65) 


By substituting, using the symbols from figure 16 equation (65) becomes: 
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Me = Me, ) (he(2) - (he(3) + Wo(3))/(h(2) - he(2)) (equation 66) 
leg = (1 - Mey - Me.) (hg(3) - h—(4) + Wo(4))/(h(3) - he(3)) (eq. 67) 
ieee (1 - My, - Me - Mag) (ng (4) - (he(S) + Wh(5))/(n(4) - he(4)) (eq. 68) 


landm=l-m-m -m. 
e 27 eo e, 


The term Wy which appears in these three equations is the enthalpy 


rise associated with the pump work. The subscript number indicates 


where m 


the temperature level at the pump inlet. The pump work, assuming the 
specific volume change across the pump is negligible is: 

= ~ 144/778 BIU/1b j 
W, (n) ve(n) Ce P )/Np X / /\bm (equation 69) 
where n is the temperature level at the pump inlet, and np is the pump 


isentropic efficiency. 
The heat added to the ammonia cycle is the product of the mass 
flow rate and the enthalpy rise in the ammonia boiler. The heat added 


Qn NH - ny (h(1) - (h¢e(2) + Wy(2)) (equation 70) 


to the ammonia cycle is also the heat rejected by the steam cycle. 


Combining equations (54) and (70), 


QA NH3 = Q% steam(l-nth steam) = m(h(1) = (he(2) + Wp(2)). (eq. 71) 
The left hand side of this equation is fixed by the steam condensing 
temperature To. The enthalpies on the right side are functions of the 
ammonia boiling and condensing temperatures, TF(1) and TF(5) respectively. 
Therefore, once these three variables are fixed, m; the mass flow rate 
through the ammonia boiler is also fixed. Recall that my is an important 
consideration in the heat transfer analysis of the steam condenser/ammonia 


boiler. 


spo) 





The thermal efficiency of the ammonia cycle is: 


4 5 
2 - 5 m: W 
n = = —S————__— 
th NH3 —— CCtt—~* Q 
A NH A NH 
c. The Combined Cycle 
The thermal efficiency of the combined cycle, kee is the total 
of the net power outputs for both cycles divided by the heat added 


to the cycle: 


+ P 
"th cc = Pst NH3 (equation 73) 
Qn steam 

recall that Pot = Qn cteam X "th Res a (equation 74) 

and P - X equation 75 
NH3 = [A NH3 % "th NH, (eq 

where Qa NH2 = QM steam (1 - nth ee, (equation 76) 

therefore, nth cc * "steam * (1 - nsteam)Nth NH3 (equation 77) 


There are several other factors which cause the overall efficiency 
to be lower than the thermal efficiency. The first is the pumping 
requirements for the ammonia condenser discussed in the preceding 
chapter, and labelled PP. The second consideration is the power 
loss associated with converting the mechanical power of the turbines 


to electrical power, and may be expressed as n_, the generator efficiency. 


g* 
The third factor is the power that must be expended to prepare the coal 
and to drive the auxiliary equipment in the plant. This factor is 


expressed as ng, the auxiliary efficiency. 
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The fourth and final factor is the boiler efficiency np which is 
a measure of how much of the energy released in the combustion of 
the fuel actually is transferred to the working fluid. 

The cycle efficiency Ne» is the efficiency of the cycle including 


the condenser pumping power: 


m 2 
Me = (Pst + Pung) - PP = np. - PPY QM steam (equation 78) 
Qa steam 


Thus the overall plant efficiency is: 
Noq = Nc X Ng X Ny X Np (equation 79) 
In this chapter, the thermodynamic analysis of the combined cycle 


was undertaken. In the next chapter, the thermodynanic and component 


considerations are combined in the economic analysis. 
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VI. ECONOMIC ANALYSIS 

The economic factors of concern in the analysis are the 
Capital and operating expenses. For simp!icity, only the cost 
differences between the combined cycle plant and ‘Bull Run' 
will be considered. 

The significant capital expenses are: 


(1) Turbomachinery - specifically the ammonia and LP 
Steam turbines. 


(2) Heat exchangers - specifically the steam condenser/ 
ammonia boiler. 


(3) Financing costs 

The significant operating expense is the fuel cost which is 
derived from the cycle overall efficiency. 

The major difference between 'Bull Run' and the combined cycle 
in turbomachinery costs jis the addition of the ammonia turbine and 
the removal of the jast stage or two of the low pressure steam 
turbine. 

Due to the large specific volume of steam at condensing conditions 
in conventional power plants, the last few stages of the LP steam 
turbine are quite large and costly. Data from DeLaval !® indicates 
that the cost of a steam turbine varies approximately linearly with 
power output as the backpressure is increased. This data fixes the 
cost variation with power at $31 per kilowatt of mechanical power output. 


The vapor pressure of ammonia at condensing conditions is much 
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higher than that of steam. Consequently, the specific volume of 
ammonia is 100 - 300 times lower than steam, causing the ammonia 
turbine to be smaller and cheaper than the steam turbine it replaces. 
The selection of salen: for the ammonia turbine is consistent with 
previous studies and existing NH3 turbine costs corrected to 1976 
dollars. 

The net power output of the combined cycle under consideration 
is the same as that of the ‘Bull Run' plant scaled up to 1000 MW. 
Therefore, any power lost by the steam turbine by the addition of the 
bottoming cycle must be compensated for by the ammonia turbine. 

The difference in turbine costs between the binary-cycle plant 

and ‘Bull Run’ may be expressed as the product of the NH3 turbine 

power and the difference in cost per kilowatt between the steam and 
ammonia turbines. Thus, 

TURBOMACHINERY COST SAVINGS = PNH3 x ($31/kw - $5/kw), (equation 80) 
where PNH3 was derived as a function of the operating parameters in 

the last chapter. 

The most significant heat exchanger cost in the combined cycle is 
that of the steam condenser/ammonia boiler. The total area of the SC/AB 
is the sum of the boiling region area, Ap sand the non-boiling area 
Aug, introduced in chapter four. 


AT = Ag + Anp (equation 81) 


Utilizing heat exchanger cost information@2, the cost of the SC/AB 
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per unit area is assumed to be $8/ft*. The cost of the steam condenser/ 
ammonia boiler is: 
COST SC/AB = A, x $a/ft® (equation 82) 
Combining equations (80) and (82) the additional initial capital 
expense of the combined cycle is: 
INITIAL CAPITAL COST = Ay x $8/ft° - Pyyg x $26/kw (equation 83) 
In practice, the funds for capital equipment are borrowed, 
incurring finance charges. 
For a loan where the annual repayment rate is the same for every 


year, the total cost of capital equipment including finance charges is: 


TOTAL COST _ INITIAL CAPITAL . i (1 + i)" . 
OF EQUIPMENT COST x Wea xn (equation 84) 


where n is the loan duration and i is the interest rate per annum. 

It must be noted that the turbomachinery cost per kw and heat 
exchanger cost per square foot of surface area represent 1976 values. 
The analysis and optimization procedure is easily modified to reflect 
price changes. 

The operating costs for the combined cycle when compared with 
those of the unmodified 'Bull Run' plant are reflected in the differences 
in fuel consumption between the two cases. The heat input to a cycle 
is the product of the fuel mass flow rate and the heating value of the 


fuel: HEAT INPUT = Mega] X HAVE a] (equation 85) 


The power output is the product of the heat input and the overall 


efficiency, POWER = HEAT INPUT/no, (equation 86) 
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Combining equations (85) and (86), the coal mass flow rate is: 


_ POWER OUTPUT (equation 87) 


m coal 
Noa X HHVcoal 


The difference in coal flow between the ‘Bull Run’ plant and 


the combined cycle is: 


_ (POWER OUTPUT)pe _ (POWER OUTPUT).. (equation a8 


MBULL RUN ~ Mec Te 
NoA X HAVcgal Non  * HAVcoal 
BULL 
RUN 








CC 


Since the analysis has assumed that the steam portion of the 
combined cycle except for the condensing region is identical to 
‘Bull Run' and that both plants have the same net power, equation (88) 
may be rewritten as follows: 
Ma p 7 Moe = Geek l (equation 89) 
coal X Na X Ng X Nb orn (me) e, 


where n,, Ng» and n, refer to equation (79) and n, is the cycle 


g 
efficiency, a modification of thermal efficiency which includes condenser 
pumping requirements. 

Equation (89) indicates that if the combined cycle efficiency 
(n.) 


cc 1S higher than the ‘Bull Run’ cycle efficiency (nc )pp fuel will 


be saved. 


The economic gain of reduced fuel consumption is: 


Ibm (eq. 90) 


DOLLAR SAVINGS PER HOUR = COAL COST $/1bm x (m,, - Mc) Ar 


BR 


Equation (90) could be used to extrapolate the savings over the plant 
lifetime, however the term 'power output' appears in equation (89). 


It would be too optimistic to design a combined cycle based on 100% 
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load demand over the entire plant life. The quantity load factor, 
LF, expresses the ratio of the expected average power output of the 


cycle to its maximum capacity. Equation (89) becomes: 


mgr - mn~.. = MAXIMUM POWER OUTPUT x LF i - 1 (eq. 91) 
EC coal xX Na xX Ng xX Nb cpp lees 





Then the fuel savings in dollars over the plant lifetime is: 


LIFETIME 2 — (eqs, 92) 
FUEL SAVINGS COAL COST $/Ton x (Mpp = "cone or Ton : x 8760 slit n yrs 


where n is the expected plant lifetime in years, as well as the loan 
duration, 

Combining the operating and capital expenses into one equation, 
the life cycle benefit in dollars of the combined cycle over the 


unmodified ‘Bull Run' plant is: 


NET = _ LIFETIME - TOTAL COST 
SAVINGS = FUEL SAVINGS — OF EQUIPMENT (FROM EQUATION 84) 


(eq. 93) 

The purpose of the optimization procedure is to maximize the 
'NET SAVINGS' of equation (93), and determine the operating parameters, 
in terms of the three independent variables Ts, TF(1) and TF(5), which 


yield the maximum economic benefit. 
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VII. OPTIMIZATION PROCEDURE 

In this chapter, the procedure for the selection of the 
optimum design bottoming cycle based on economic and engineering 
factors will be outlined. 

The optimum design may be specified in terms of three operating 
parameters previously identified as the independent variables in the _ 
analysis, they are: 

(1) The steam condensing temperature, Ts. 

(2) The ammonia boiling temperature, TF(1). 

(3) The ammonia condensing temperature, TF(5). 

Once values for these variables are chosen, a thermodynamic 
analysis of the combined cycle is undertaken which yields the combined 
cycle thermal efficiency and the power split as the output. Power 
Split is defined as the amount of power produced by each of the two 
cycles. ’ 

Next, the heat transfer analysis for the steam condenser/ammonia* 
boiler specifies the size of the heat exchanger. The ammonia Paaeneee 
is then analyzed with the condenser pumping power as the result. 

The information derived from these analyses is combined in 
equation (93), which combines capital and operating costs into a value of 
net savings over the plant lifetime. 

A flow diagram of the optimization procedure is shown as figure 


17. Note that the value of net savings resulting from the analysis is 
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very sensitive to the economic factors used in the cost computations. 
The most significant factors are fuel and heat exchanger costs. 

Thus far, the optimization has produced a value of net savings 
of the combined cycle plant over that of the ‘Bull Run' plant scaled 
up to 1000 MW for one set of the three variables Ty, TF(1), and TF(5). 

In order to determine which set will yield the optimum plant, the entire 
analySis must be carried out for all possible combinations of the three 
independent parameters within the range of values where the bottoming 
cycle is expected to operate. 

The analysis may be accomplished independent of environmental 
conditions with regard to selection of values of the steam condensing 
temperature, TS, and the ammonia boiling temperature TF(1). However, 
the selection of the ammonia condensing temperature, TF(5) is very much 
affected by the temperature of the condenser coolant water at the inlet. 

As discussed in chapter four, the relationship between the condenser 
water inlet temperature and the ammonia condensing temperature is restricted by 
considerations of pumping power and tube corrosion resulting from high 
water mass flow rates and velocities. 

Utilizing equations (13) and (78) the cycle efficiency may be 
calculated for various values of temperature difference between the con- 
densing ammonia and the inlet water for fixed parameters Ts, TF(1), 
and TF(5). Additionally, the water velocity through the tubes as a func- 
tion of the same parameters may be calculated. The results of these 
calculations are plotted as figure 18. The values of steam condensing 


temperature, ammonia boiling temperatures, and ammonia condensing tem- 
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FIGURE 18 CONDENSER WATER VELOCITY AND CYCLE EFFICIENCY VS. DIFFERENCE 
BETWEEN NH, CONDENSING TEMPERATURE AND WATER INLET TEMPERATURE 
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perature used for figure 18 are representative, near optimum values. 
Notice that the cycle efficiency is at a maximum at approximately a 

15°F temperature difference between the condensing ammonia and the 

inlet coolant water. Furthermore, the same temperature difference 
yields a water velocity of 12.5 ft/sec, close to the corrosion limit 

of 15.0 ft/sec. Thus it appears that the optimum selection of TF(5), 
the ammonia condensing temperature is 15°F above the available water 
source temperature derived from environmental information. This conclu- 
sion agrees with earlier condenser design studies .<' 

Replacing the ammonia condensing temperature TF(5) as one of the 
independent parameters by the water inlet temperature Tin and keeping 
the 15°F difference between TF(5) and Tyy for all possible designs, the 
optimization computations for different sets of the independent par- 
ameters may now proceed. 

Since the combined cycle under consideration is compared to the 
‘Bull Run' plant, the environmental conditions of the ‘Bull Run' plant 
are used. At the ‘Bull Run’ site, the average condenser inlet temper- 
ature is cere and this value is incorporated into the optimum design 
analysis. However, other values of condenser inlet temperature are 
included in the results. 

Therefore, only the steam condensing and ammonia boiling tem- 
peratures need be varied to find the optimum design. In practice, 
the steam condensing temperature is fixed and the temperature difference 
between the condensing steam and boiling ammonia is varied. To present 


the results on a non-dimensional basis, this temperature difference is 


=712 





expressed as the difference between the steam condensing temperature 
and the entropy averaged temperature of heat addition to the ammonia, 
T, divided by the steam condensing temperature, 


ete (1. - 7)/ T. (equation 94) 


S 
where T is a function of the ammonia boiler inlet and outlet conditions. 

In order to carry out the optimization procedure, the many 
engineering and economic factors previously discussed must be given 
values. The optimization computer program (see appendix one), will 
carry out the analysis for any values of these factors. This flexi- 
bility enables this optimization to be used for future bottoming cycle 
design utilizing up-to-date economic considerations. The symbols and 
values used in this analysis are shown in table 4. 

Thus far, the derivation of the analysis and the procedure for 
identifying the optimum design combined cycle have been presented. 
In the next chapter, the results of the optimization calculations are 
presented along with the specifications of the optimum bottoming cycle 


compatable with the ‘Bull Run’ plant. 
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EFFICIENCIES 
Auxiliary Efficiency 
Boiler Efficiency 
Bull Run Cycle Efficiency 
Generator Efficiency 

Pump Efficiency 
Turbine Isentropic Efficiency 
Turbine Mechanical Efficiency 
(Includes exit, stage and 
mechanical losses) 

HEAT EXCHANGER SPECIFICATIONS 


Ammonia Boiler Tube Diameter 


Ammonia Boiler Tubes in Vertical 
Bank 


Ammonia Condenser Tube Diameter 


Ammonia Condenser Tubes in Vertical] 
Bank 


ECONOMIC FACTORS 
Coal Cost (111. #6) 
Heat Exchanger Cost 
Higher Heating Value (I11. #6) 
Interest Rate 
Load Factor 
Loan Duration and Plant Life 
Turbine Cost (Ammonia) 


Turbine Cost (LP Steam) 


TEXT 
SYMBOL 


Me 


HHV 


EE 


COMPUTER 
SYMBOL 


NA 
R4 
R5 
R7 
NP 
NS 
NM 


D2 
NOT 


D3 


INCLUDED IN 
A CONSTANT 


R6 
RI 


IR 
R8 
R9 
R3 
R2 


TABLE 4 SUMMARY OF ENGINEERING AND ECONOMIC FACTORS 
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VALUE 


94% 
89% 
45.78% 
She ey 
90 &% 
84 = &% 
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son0" 


200 


$50/Ton 
$ 8/ft* 


10788. BIU/1bm 


12% ; 

8 

30 years 
$ 5/kw 


$31/kw 





VIII. OPTIMIZATION RESULTS 

Utilizing the optimization procedure outlined in the last 
chapter, the thermal efficiency and lifetime net savings for each 
set of combined cycle parameters are calculated. Assuming the 
average condenser water inlet temperature to be 55°F, and 
incorporating the engineering and economic factors enumerated in 
table 4, an optimum design may be predicted. 

Figure 19 is a plot of the combined cycle efficiency versus 
the non-dimensional temperature difference AT/Ts for several different 
values of Ts, the steam condensing temperature. As expected, the cycle 
efficiency increases as the temperature difference decreases. 
Depending on the value of AT/Ts, either a steam condensing temperature 
of 131°F or 147°F yields the maximum efficiency. This efficiency peak 
is caused by wetness considerations discussed in chapter 4. 

However, as the non-dimensional temperature difference decreases, 
the area of the steam condenser/ammonia boiler increases. The 
opposing trends of efficiency and heat exchanger size variation with 
the non-dimensional temperature difference are the most significant factors 
in the optimization routine. Figure 20 shows the steam condenser/ 
ammonia boiler areas versus the non-dimensional temperature difference 
for several values of Tos the steam condensing temperature. Note that 
the heat exchanger area is slightly smaller for T. = 147°F than for 


131°F. This trend is caused by the change of the steam properties which 
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FIGURE 19 COMBINED CYCLE EFFICIENCY VERSUS NON-DIMENSIONAL 
TEMPERATURE DIFFERENCE FOR VARIOUS VALUES OF 
STEAM CONDENSING TEMPERATURE. WATER INLET 
TEMPERATURE = 55°F. 
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FIGURE 20 STEAM CONDENSER/ AMMONIA BOILER HEAT TRANSFER AREA 
VERSUS NON-DIMENSTONAL TEMPERATURE DIFFERENCE FOR 
VARIOUS VALUES OF STEAM CONDENSING TEMPERATURE. 
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affect the steam side heat transfer coefficient (seeequation (23)) as 
the condensing temperature varies. 

Combining the effects of cycle efficiency, heat exchanger 
size, and the other economic factors in the cost equation (equation (93)), 
a plot of lifetime cycle savings is presented as figure 21. 

It is obvious that the maximum savings occurs where the steam 
condensing temperature Tc = 147°F, and the non-dimensional temperature 
difference AT / Ts = .09. This value of AT / Ts corresponds to an 
ammonia boiling temperature TF(1) = 142°F. 

The predicted lifetime combined cycle savings over the ‘Bull Run' 
plant is $98 million or about .047 cents/kw - hr. The combined cycle 
efficiency is 47.2%, about 1.4% higher than the ‘Bull Run' plant 
with an initial additional capital cost of only $9.1 million. 

The flow chart and thermodynamic specifications of the combined 
cycle are presented as figure 22 and table 5 respectively. 

In order to determine the combined cycle performance in the off- 
design condition, several assumptions must be made. They are: 


(1) Maintain steam condensing and ammonia boiling 
temperature at design points. 


(2) Maintain a 15°F temperature difference between the 
ammonia condensing temperature and the coolant water 
inlet temperature. 


(3) Utilize DeLaval estimate of off-design ammonia turbine 
performance (figure 23). 


Carrying out the off-design analysis, incorporating the exit 
loss criteria from figure 23, a plot of combined cycle efficiency versus 


coolant water inlet temperature is shown as figure 24. For water inlet 
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FIGURE 21 PLANT LIFETIME NET SAVINGS VERSUS NON-DIMENSIONAL 
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Fue | 
Fue] 


Air 
Exhaust 


Ash 
H90 


H50 
H50 
Ho0 
H90 


Ho0 


NH3 
NH3 
NH3 
NH3 


H50 
H50 


NO. 


7 

Qk 

9 
10 
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DESCRIPTION 


Raw fuel (13% H0) 
Prepared fuel 

(6% Ho0) 
Combustor inlet air 


Products of 
combustion 


Ash 


Water extracted 
from fuel 


Condenser exit 
Boiler inlet 


Boiler exit/HP 
turbine inlet 


HP turbine exit/ 
reheater inlet 


Reheater exit/IP 
turbine inlet 


IP turbine exit/ 
condenser inlet 


NH3 condenser exit 
NH3 boiler inlet 


NH3 boiler exit 
NH3 turbine exit 


Coolant water inlet 


Coolant water exit 


* INCLUDES FEEDHEATING 


TABLE 5(PART 1) 
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THERMODYNAMIC SPECIFICATIONS OF THE COMBINED CYCLE 
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COMPONENT POWER HEAT TRANSFER 


Steam boiler 1876.0 MW 
Steam reheater 400.5 MW 
HP steam turbine 272.2 MW 
IP steam turbine 682.6 MW 
Steam condenser/NH3 boiler 1321.7 MM 
Ammonia turbine 123.9 MW 
Ammonia condenser 1201.9 MW 


Gross total power 1078.7 MW 


Generator losses -10.8 
Auxiliary power ~64.1 
Condenser pumps - 3.8 
NET POWER 1000.0 MW 


Combined cycle: Thermal efficiency = 47.38% 
Cycle efficiency = 47.2% 
Overall efficiency = 39.1% 


Steam cycle : nep = 41.94% 
Ammonia cycle : NNH3" 9.37% 


TABLE 5(PART 2) THERMODYNAMIC SPECIFICATIONS OF THE COMBINED CYCLE 
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FIGURE 24 OFF-DESIGN PERFORMANCE OF COMBINED CYCLE 
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temperature above 55°F, the pressure ratio used in figure 23 becomes 
greater than one. For pressure ratios greater than unity, the exit 
loss is assumed to be 1%. 

Optimization results for condenser water inlet temperatures of 
35°F, 45°F, and 65°F are shown in figure 25 through 30. The parameters 
of the optimum design for each value of condenser inlet temperature are 


shown in table 6. 


CONDENSER STEAM NH3 NH 
INLET CONDENSING CONDENSING — BOILING NET 
TEMPERATURE TEMPERATURE TEMPERATURE TEMPERATURE SAVINGS 
35°F 131°F 50°F 126°F $233 million 
45°F 147°F 60°F 142°F $162. million 
55°F 147°F 70°F 142°F = $98 million 
65°F 147°F 80°F 142°F = $- 29 million 


TABLE 6 OPTUMUM DESIGN PARAMETERS 


The net savings for each optimum design taken from table 6 
is plotted in figure 31. As expected, the combined cycle is most 
attractive where the ambient temperature is lowest. More significantly, 
figure 3] indicates at what value of condenser inlet temperature the 
addition of a bottoming cycle is no longer economically viable, approximately 
69°F. 

In this chapter, information has been presented to permit the 
Selection of an optimum ammonia bottoming cycle with the available water 


Source temperature as the only input variable. 
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FIGURE 25 COMBINED CYCLE EFFICIENCY VERSUS NON-DIMENSTONAL 
TEMPERATURE DIFFERENCE 
CONDENSER INLET TEMPERATURE = 35°F 
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FIGURE 26 COMBINED CYCLE NET SAVINGS VERSUS NON-DIMENSIONAL 
TEMPERATURE DIFFERENCE 
CONDENSER INLET TEMPERATURE = 35°F 
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FIGURE 27 COMBINED CYCLE EFFICIENCY VERSUS NON-DIMENSIONAL 
TEMPERATURE DIFFERENCE 
CONDENSER INLET TEMPERATURE = 45°F 
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FIGURE 28 COMBINED CYCLE NET SAVINGS VERSUS NON-DIMENSTIONAL 
TEMPERATURE DIFFERENCE 
CONDENSER INLET TEMPERATURE = 45°F 
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FIGURE 29 COMBINED CYCLE EFFICIENCY VERSUS NON-DIMENSIONAL 
TEMPERATURE DIFFERENCE 
CONDENSER INLET TEMPERATURE = 65°F 
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IX. SUMMARY AND CONCLUSION 

a. SUMMARY 

The purpose of a bottoming cycle is to improve the efficiency 
of a power plant by enabling the temperature of heat rejection to be 
lower than what is practical for a steam plant. 

By improving the cycle efficiency, a bottoming cycle equipped 
power plant will consume less fuel, reject less heat to the environment, 
and deposit fewer harmful products of combustion into the air. 

Ammonia seems to be the best choice as the working fluid for a 
bottoming cycle designed to operate in sub-position to a conventional 
steam plant such as ‘Bull Run’. 

The bottoming cycle is defined by three operating parameters: 

(1) Ts, the steam condensing temperature. 

(2) TF(1), the ammonia boiling temperature. 

(3) TF(5), the ammonia condensing temperature. 

In order to determine the operating parameters of the optimum 
bottoming cycle, a thermodynamic analysis was conducted to find the cycle 
efficiency as a function of the operating parameters. Next, the steam 
condenser/ammonia boiler and turbomachinery sizes were calculated. 

These operating and capital expenses are combined into an equation which 
yields the net lifetime savings of a combined cycle plant over the 
unmodified steam plant. The set of operating parameters which gives the 


maximum economic benefit is taken as the optimum design. 
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The steam condensing temperature which gives maximum efficiency 
is that value which minimizes the total power lost due to wetness in 
both the steam and ammonia cycles. 

Limitations on the ammonia condensing temperature are due to cor- 
rosion and pumping power considerations on the water side of the condenser. 
The ammonia condensing temperature for optimum performance is 15°F 
above the condenser water inlet temperature. 

Thus, the economic analysis hinges on the selection of the ammonia 
boiling temperature, or more precisely the difference between the 
steam condensing temperature and the entropy averaged temperature of 
heat addition to the ammonia. As this temperature difference decreases, 
the thermodynamic efficiency of the combined cycle increases but the 
area of the steam condenser/ammonia boiler increases. The economic 
analysis indicates what value of this temperature difference makes the 
addition of an ammonia bottoming cycle most attractive. 

Given the present economic factors enumerated in table 4, the 
optimum bottoming cycle to be installed at the ‘Bull Run' site has 


the following operating and economic parameters: 


(1) Steam condensing temperature Ts = 147°F 
(2) Ammonia boiling temperature TF(1) = 142°F 
(3) Ammonia condensing temperature TF(5) = 70°F 
(4) Condenser water inlet temperature Tjy = SS°F 
(5) Cycle net power p = 1000 MW 
(6) Overall cycle efficiency Non = 39.1% 


ZOa= 





(7) Cycle efficiency increase over An = lean 


‘Bull Run' 
(8) Lifetime net savings $98 million or .047¢/kw-hr 
(9) Added initial capital cost $9.1 million 


b. CONCLUSION 

The ammonia bottoming cycle is economically and technologically 
feasible using existing technology and based on current economic 
conditions, provided the condenser water inlet temperature is less 
than 69°F. If the cost of fuel continues to increase, which is likely, 
the bottoming cycle becomes even more attractive. 

The addition of a bottoming cycle seems to be suitable for large 
conventional utility steam plants, and pressurized-water type nuclear 
plants employing water-cooled condensers. However, due to the large 
Size of the steam condenser/ammonia boiler which is on the order of 
1.5 x 106 Square feet of surface area, the bottoming cycle is not deemed 
practical for applications where power plant size is a major consideration. 
Additionally, the possibility of a major ammonia leak precludes the 


installation in a closed-in environment such as a ship. 
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APPENDIX I THE COMPUTER PROGRAM 
a. INTRODUCTION - The computer program used in the optimization/ 
analysis was run on the M70 and M80 computers at the joint 
Mechanical/Civil Engineering computer facility at MIT. The programming 
language used was Fortran IV. 
b. IDENTIFICATION OF COMPUTER VARIABLES 
VARTABLE REPRESENTATION 


A Intermediate variable which sums AB's 
AB Area of one quality region of boiling section in SC/AB 
ANB Area of non-boiling section of SC/AB 
AT Total area of steam condenser/NH3 boiler 
A2 NH3 condenser area 
A5 Difference between density of saturated liquid and 
Saturated vapor of NH3 at condensing temperature 
A6é Density of saturated liquid NH3 at condensing temperature 
A7 Difference between NH3 condensing temperature and the 
tube wall temperature 
BETA Part of condensing heat transfer coefficient 
BI Average quality in a section of the boiling region of 
NH3 boiler 
c Heat transfer coefficient on boiling side of SC/AB 
(first approximation) 
CAPT Additional capital cost of bottoming cycle 
CP specific heat at constant pressure of NH3 
C1 Heat transfer coefficient on boiling side of SC/AB( final value) 
C3 Heat transfer coefficient on NH3 side of non-boiling section 
of SC/AB 
C4 Heat transfer coefficient on steam side of non-boiling 
section of SC/AB (first approximation) 
C5 Error between C3 and C4 in iteration loop 
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VARIABLE 


C6 
D 


DELTA 


DP 
DI 


D2 
D3 


EFF 
EFFOA 


EFFST 


ENTH 
ETA 
E] 


E2 
ES 


FI( ) 
it) 


GW 


H( ) 
HC 


HCOND 
HF( ) 
HFG( ) 


REPRESENTAT ION 
Final value of steam condensing heat transfer coefficient 


Heat transfer coefficient on condensing side of SC/AB 
(first approximation) 


Difference between the steam condensing temperation (TS) 
and the average temperature of NH3 heat addition (TBAR) 


Pressure drop on water side of NH3 condenser 


Heat transfer coefficient on condensing side of SC/AB 
(final value) 


NH3 boiler tube diameter 

NH3 condenser tube diameter 

Error in heat flux calculation in SC/AB (C-D) 
Thermal efficiency of NH3 cycle 


Thermal efficiency of combined cycle excluding condensing 
pump requirements 


Efficiency of steam cycle at specified steam condensing 
temperature 


Corrected enthalpy used in NH3 condensing calculation 
Overall thermal efficiency of combined cycle 


Absolute difference between first and second approximation 
in NH3 condenser water outlet temperature routine 


Error in first approximation of heat flux in NH3 condenser 
Absolute value of E2 

Chen factor for boiling 

Intermediate function used to calculate HG( ) 
Intermediate function used to calculate SG( ) 

NH3 mass flow per unit area in NH3 boiler 

Mass flow rate/unit area of water in NH3 condenser 
Enthalpy (NH3) 


Convective component of the heat transfer coefficient in 
the NH boiler 


Condensing NH3 heat transfer coefficient 
Enthalpy of saturated liquid (NH; ) 
Latent heat of evaporation (NH3) 
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VARIABLE 
HG( ) 


HNUB 


HS( ) 


HW 


IR 


KF 
KF] 


KW 


LMTD 


ME( ) 


MNH- 
MW 


NM 


NOT 
NP 

NS 

P( ) 
PF 

PG 
PH20 
PNH3 
PP 
PUMPW 
PX( ) 


QA 


REPRESENTATION 

Latent heat of evaporation vapor (NH, ) 

Nucleate boiling component of heat transfer coefficient 
in NH, boiler 

Enthalpy (based on isentropic expansion in NH turbine) 


Convective heat transfer coefficient of water side 
of NH. condenser 


Counting integer in DO LOOP 
Capitalization rate 


Thermal conductivity of NH, in NH3 boiler 
Thermal conductivity of NH, in NH, condenser 
Thermal conductivity of water 

Temperature step side in AT variation routine 
Log mean temperature difference SC/AB 

NH. turbine mass flow extraction for feedheating 
Mass flow rate through NH3 boiler 

Mass flow rate through NH, condenser 

Counting integer in DO LOOP 


Turbine efficiency after mechanical, exit and stage 
losses 


Number of tubes in NH, boiler 

Pump efficiency 

Isentropic efficiency of turbines 

Pressure (in thermodynamic properties of NH3 sub routine) 
Density of NH3 Saturated liquid 

Density of NH, Saturated vapor 

Power produced by steam cycle 

Power produced by NH, 

NH3 condenser water pumping power 

Total NH, feed pumping power 

Intermediate function used to calculate P( ) 
Heat added to NH cycle 


Enthalpy added to NH, in SC/AB 
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VARIABLE 


QAD 
QB 
QNB 
QR 
RETP 
RI 
R2 


R3 
R4 


R5 
R6 
R7 
R8 
R9 

S( ) 

SAVE 


SE ) 
SG(_ ) 
ST 

S] 

T( ) 
TAV 
TBAR 
TF( ) 
TIN 
TOUT 


TOUTB 


TS 
TURBW 


REPRESENTATION 

Heat added to steam cycle 

Total heat added to boiling section of NH, boiler 
Total heat added to non-boiling section of NH. boiler 
Total heat rejected by combined cycle 

Two-phase flow Reynolds number 

Heat exchanger cost $/FT 

Steam turbine cost $/MW 

NH. turbine cost $/MW 

Boiler efficiency 

Bull Run thermal efficiency 

Coal cost $/TON 

Generator efficiency 

Load factor 

Plant life and loan duration 

Entropy 


Lifetime savings gained by bottoming cycle equipped 
plant over base case 


Entropy of saturated liquid (NH3) 

Entropy of saturated vapor (NH) 

Surface tension of NH3 in NH boiler 

Boiling suppression factor in two-phase flow 
Temperature (°R) 

Average water temperature in NH, condenser 
Average temperature of heat addition of NH, 
Temperature (°F) 

Inlet water temperature in NH, condenser 


Outlet water temperature in NH3 condenser (first approx- 
imation ) 


Outlet water temperature in NH, condenser (second approx- 
imation ) 


Steam condensing temperature (°F) 


Total hale turbine power 
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VARIABLE 
TT 
TW 
U 
UB 


UC 
UF 
UF] 
UG 
UNB 


UW 
VF( ) 
VFG( ) 
Vet ) 

VW 
WP( ) 
WT( ) 

X( ) 
XS(_) 

XTT 

YTT 


REPRESENTATION 

NH. condenser tube wall temperature 
SC/AB tube wall temperature 

Overall heat transfer coefficient 


Overall average heat transfer coefficient in boiling 
section of the SC/AB 


Overall heat transfer coefficient in NH condenser 

Absolute viscosity of NH3 Saturated liquid in NH, boiler 
Absolute viscosity of NH, saturated liquid in NH, condenser 
Absolute viscosity of NH, Saturated vapor in NH, boiler 


Overall heat transfer coefficient in non-boiling section 
of the SC/AB 


Average absolute viscosity of water in NH3 condenser 
Specific volume of NH_ saturated liquid 

VG{ ) - VF( ) 

Specific volume of NH. Saturated vapor 

Velocity of water in NH3 condenser 

NH feed stage pump work 

NH turbine stage work 

Quality 

Quality based on isentropic expansion in NH3 turbine 
Martinelli parameter for boiling 

Inverse of XTT 

% Quality in NH3 boiler section (Bl x 100) 
Delta/Ts(°R) 
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c. PROGRAM DESCRIPTION 

For clarity, the computer program on the following pages is 
divided into nine zones, labelled 'A' through 'I'. Each zone 
represents a distinct function in the optimization routine. 

Zone 'A' includes the identification of variables as either 
integers or real numbers. The ‘dimension’ statements reserve core 
space for the subscripted variables. 

Zone 'B' consists of the engineering and economic factors 
listed in table four. It is these factors which the economic 
analysis hinges upon. Accordingly, these parameters are easily 
changed to allow an up-to-date analysis. Also included in zone 
'B' are the ammonia condensing temperature and the coolant water 
inlet temperature for the plant under consideration. 

The steam condensing temperature, steam cycle efficiency, 
and part of the steam side heat transfer coefficient are included 
in zone 'C'. Also included are the temperature difference between 
the condensing steam and boiling ammonia, as well as the pump and 
turbine efficiencies. 

Part 'D' establishes the three intermediate temperature levels 
in the ammonia cycle and generates the thermodynamic properties of 
ammonia needed for the analysis. The equation of state for the ammonia 


are empirical curves, fitted to the tabulated data’’. 
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The thermodynamic analysis of the ammonia cycle and the 
calculation of the combined cycle thermal efficiency occurs 
in zone 'E'. It is in this section where the non-dimensional 
temperature difference between the condensing steam and boiling 
ammonia is calculated. 

In zone 'F' the heat transfer properties of the ammonia are 
calculated from empirical Pe arion ee 

Zone 'G' is the most complex part of the entire program. In 
this section the steam condenser/ammonia boiler area is calculated. 
Note that in statements 55 and 77 the tube wall temperature is 
incrementally changed. As stated in chapter four, the heat flux 
equation for both sides of the steam condenser/ammonia boiler can only 
be solved by interating until the wall temperature which satisfies 
the equality is found. In order to minimize the running time of the 
program, an exact equality of both sides of equation (34) is not 
demanded, but rather it is required that the two results are within 
1% of each other. Without this relaxation of the solution, the 
temperature increment required to obtain a converging iterative solution 
would increase the running time and cost of the program by several orders 
of magnitude. Note also that the heat exchanger analysis is repeated 
for the 10 quality regions discussed in chapter four, as well as for the 
non-boiling region. 

The ammonia heat transfer properties and condenser pumping power 
are computed in zone 'H'. Once again the tube wall temperature is found 


by iteration. 
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In zone 'I', the condenser pumping power is incorporated into 
the cycle efficiency. The initial capital, operating, and lifetime 
costs are also calculated in this region. The output statements com- 


prise the final section of the program. 
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APPENDIX II STEAM CYCLE DATA 
The computer program discussed in Appendix I requires the 
following information about the steam cycle as inputs: 


(1) The steam condensing temperature, Tc. 

(2) The steam cycle thermal efficiency at Ts. 

(3) The portion of the steam side heat transfer (condensing) 

coefficient in the steam condenser/ammonia boiler, 
which is a function of T, i.e. the fluid properties. 

Carrying out the modification of the ‘Bull Run’ plant detailed 
in chapter five, the steam cycle thermal efficiency versus the steam 
condensing temperature is calculated. 

The properties involved in the calculation of the steam side 
heat transfer coefficient and their relationship are shown in 
equation (23). The portion of equation (23) which includes only the 
fluid properties is called 'B,'. 


Table A2-1 shows values of thermal efficiency, nth and 'Bo' 


steam 
for the different values of steam condensing temperature used in the 


analysis. 

*s(°F) Mths team 8,(BTU/(hr-ft-(°F)*’°)) 
1.7 "4578 3438.6 

110 4451 3642.4 

127 4336 3760.3 

131 4308 3791.7 

147 4194 3917.2 

162.2 4084 4085.5 

182.2 3952 4270.1 


TABLE A2-1 VARIATION OF STEAM CYCLE EFFICIENCY AND ‘8,' WITH THE STEAM 
CONDENSING TEMPERATURE 
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